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ABSTRACT 

We present a survey of H 2 jets from young protostars in the Vela-D molecular 
cloud (VMR-D), based on Spitzer-IKAC data between 3.6/xm and 8.0 /zm. Our 
search has led to the identification of 15 jets (2 new discoveries) and about 70 
well aligned knots within 1.2 deg 2 . We compare the IRAC maps with observa- 
tions of the H 2 1-0 S(l) line at 2.12 /iin, with a Spitzer-MIPS map at 24 /zm and 
70 /im, and with a map of the dust continuum emission at 1.2 mm. From such a 
comparison we find a tight association between molecular jets and dust peaks. 
The jet candidate exciting sources have been searched for in the published cat- 
alog of the Young Stellar Objects of VMR-D. In particular, we searched for all 
the sources of Class II or (preferentially) earlier which are located close to the 
jet center and aligned with it. Furthermore, the association between jet and ex- 
citing source was validated by estimating the differential extinction between the 
jet opposite lobes. We are able to find a best-candidate exciting source in all but 
two jets, for which two alternative candidates are given. Four exciting sources 
are not (or very barely) observed at wavelengths shorter than 24 /iin, suggesting 
they are very young protostars. Three of them are also associated with the most 
compact jets (projected length < 0.1 pc). 

The exciting source Spectral Energy Distributions (SEDs) have been constructed 
and modeled by means of all the available photometric data between 1.2 /im and 
1.2 mm. From SEDs fits we derive the main source parameters, which indicate 
that most of them are low-mass protostars. 

A significant correlation is found between the projected jet length and the [24] - 
[70] color, which is consistent with an evolutionary scenario according to which 
shorter jets are associated with younger sources. A rough correlation is found 
between IRAC line cooling and exciting source bolometric luminosity, in agree- 
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ment with the previous literature. The emerging trend suggests that mass loss 
and mass accretion are tightly related phenomena and that both decrease with 
time. 

Subject headings: Stars: formation - surveys -ISM: individual (Vela Molecular Ridge) 
- ISM: clouds - ISM: jets and outflows - infrared: stars 
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1. Introduction 

Supersonic flows observed as collimated jets are frequently associated with young 
protostars of different masses and evolutionary stages. Jets represent an essential tool by 
which the driving source removes its excess angular momentum, allowing the collapse to 
go on. Their interaction with the ambient medium occurs via excitation of atoms and 
ions, whose cooling lines represent the way through which jets are seen. Such cooling lines 
spread over a wide frequency range from optical-UV up to far-IR. In particular, molecular 
hydrogen imaging of the v = 1-0 S(l) line at 2.12 /iin is extensively used for identifying 
knots of molecular emission, while the spectroscopy of the H 2 ro-vibrational lines lying in 
the near-IR range (1-5 /zm) is largely effective both to probe the molecular gas at thousands 
Kelvin and to infer the main excitation mechanism(s) (fluorescence, C-ontinous or J-ump 
shocks, see e.g. Black & van Dishoeck 1987; Kaufman & Neufeld 1996; Hollenbach & 
McKee 1989; Smith 1995). 

Statistically, jet images are powerful in correlating the properties of the young 
population of the parental cloud (e.g. age, degree of clustering, Initial Luminosity Function 
- ILF) with their common features (e.g velocities, dynamical ages, collimation and position 
angles). Such correlations are undoubtedly fundamental for investigating how the cloud 
parameters affect the onset, the properties and the morphology of the jets and how these 
latter influence the cloud itself. 

During the last two decades, hundreds of individual Young Stellar Objects (YSOs) have 
been imaged with near-infrared narrow band filters (H 2 and [Fe II]), although extensive 
surveys able to cover square degree scale star forming regions are not so numerous (for some 
exceptions see eg. Stanke et al. 1998, Davis et al. 2008, 2009). In this sense the digital 
catalog of Molecular Hydrogen Objects (MHO, Davis et al. 2010) undoubtedly represents a 
relevant tool of investigation. 
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The InfraRed Array Camera (IRAC, Fazio et al. 2004) on board of the Spitzer Space 
Telescope (Werner et al. 2004) offers a unique chance to carry out large field searches of 
both jets and their obscured exciting sources (e.g. Cyganowski et al. 2011), Indeed, IRAC 
band-passes (between 3.6 /im and 8.0 /xm) contain numerous H 2 emission lines (see Fig.l of 
Smith & Rosen 2005) and are therefore suited to find jets and to study their properties, as 
testified by an increasing number of observational works (e.g. Teixeira et al. 2008; Smith 
et al. 2006; Velusamy et al. 2007; Neufeld & Yuan 2008, Takami et al. 2011, Ybarra et al. 
2010, Ohlendorf et al. 2012, Lee et al. 2012). Such an observational capability led to the 
implementation of detailed interpretative tools (Smith & Rosen 2005; Neufeld & Yuan 2008; 
Ybarra & Lada 2009 (YL09); Takami et al. 2010) able to analyze IRAC maps of molecular 
flows. However, while the mere jet identification is a straightforward task, the derivation 
of the physical conditions from such photometric maps has to be done with some caution 
(Noriega- Crespo et al. 2004, Neufeld & Yuan 2008). Indeed the H2 emission lines inside 
the IRAC bands come from vibrational levels corresponding to very different excitation 
temperatures (from 800 K to 4000 K), and contamination from PAH, CO and atomic lines 
can be significant (Takami et al. 2010). 

In addition to allow jets identification, mid-infrared IRAC images have a major role 
in finding the embedded exciting sources that could have escaped detection at shorter 
wavelengths. Otherwise, IRAC observations allow to probe the continuum emission of 
already identified exciting sources, but in a frequency regime where most of their energy is 
likely emitted. Particularly suited are the IRAC surveys of star forming regions for which a 
census of the young stellar population has been already taken, since they offer an immediate 
correlation between the identified jets and the nature of the exciting sources. A well suited 
test-bed is the cloud D of the Vela Molecular Ridge (hereinafter VMR-D, Murphy & May 
1991), an active site of star formation in the Galactic plane showing a rich phenomenology 
associated with star formation processes, such as outflows (Wouterloot & Brand 1999; Elia 
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et al. 2007), jets (Lorenzetti et al. 2002b; Giannini et al. 2005), and evidences of clusters 
as well as isolated protostars (Massi et al. 2000, 2003, Giannini et al. 2007, Strafella et al. 
2010, hereafter SEC10). 

In the present paper we analyze the IRAC maps of VMR-D in order to : i) obtain an 
IRAC-based census of the protostellar jets in VMR-D, it) compare the jets found in the 
IRAC images with those we have already detected in the H2 1-0S(1) at 2.12 nm and validate 
(or not) the applied IRAC method; Hi) identify new exciting sources or provide the mid-IR 
photometry for those already known; iv) derive the morphology of the jets, making use also 
of data in other bands, and v) correlate the properties of the jets with those of their exciting 
sources. 

The paper has the following structure: we first briefly present the data (Sect. 2) and 
explain how jets are identified and analyzed (Sect. 3). In Sect. 4 we search for the jet 
exciting sources and present their photometry. In Sect. 5 we discuss our results, that are 
briefly summarized in Sect. 6. Finally, we give in AppendbdAja brief description of each 
detected jet. 

2. Observational data 

2.1. IRAC, MIPS and millimeter data 

We have observed VMR-D as part of the IRAG-Spitzer Cycle 3 GTO (PID 30335, 
PI: G. Fazio). The observations cover a region of about 1.2 square degree delimited in 
longitude by 264°29' < I < 263°00' and in latitude by 0°42' < b < 0°50'. We have presented 
the photometric results in SEC10, who describe in detail the observational strategy, source 
extraction and photometry. This work resulted in a catalog, which we use in the present 
paper to identify the jet exciting sources (ES). In the same catalog it is also listed the 
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MlPS-Spitzer photometry at 24 /im and 70 /im (Giannini et al. 2007, SEC10) as well the 
1.2 mm photometry (Massi et al. 2007, referred also as SIMBA photometry), which are 
used here to construct the SED of the exciting sources. 

2.2. H 2 2.12 /im images 

To compare the Spitzer maps with the H 2 emission in the near-infrared, we use images 
of the H2 1-0 S(l) ro-vibrational transition at 2.12/im we obtained on December 2001 with 
SofI (Lidmann et al. 2006) at ESO-NTT (La Silla, Chile). These images cover the peaks of 
the dust continuum emission found by Massi et al. (2007) and have been already searched 
for evidences of shocked emission by Lorenzetti et al. (2002b) and De Luca et al. (2007). 
Out of 18 H2 SofI fields, we have detected shocked emission in the 11 fields displayed in 
Figure[TJ which are overlaid on the IRAC image of VMR-D at 4.5 /im (left panel). In the 
same Figure, right panel, the dust emission image at 1.2 mm is depicted, as well. In one 
case (jet # 15) the 2.12 /mi image has been obtained with ISAAC (Cuby et al. 2004) at 
ESO-VLT (Paranal, Chile), and reported in Giannini et al. (2005). 

2.3. Complementary data 

2.3.1. Public surveys data 

In order to supplement the photometry of the ES of the jets in the wavelength ranges 
adjacent to those probed with IRAC/MIPS, we used the 2MASS (Skrutskie et al. 2006) 
photometry in the J, H and K s bands, and the recently delivered data of the Wide-field 
Infrared Survey Explorer (WISE, Wright et al. 2010) taken at 3.4 /im, 4.6 /im, 12 /im, and 
22 /im. Finally, data obtained at 250 /mi, 350 /im, and 500 /im with the Balloon-borne 
Large-Aperture Submillimeter Telescope (BLAST, Pascale et al. 2008) and compiled by 
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Netterfield et al. (2009) and Olmi et al. (2010) have also been considered in the SED's 
building. 



2.3.2. CO APEX data 

12 CO(3-2) maps of four jets (namely # 1, 2, 3, 5) were obtained from September to 
December 2006 with the APEX-2A heterodyne receiver, mounted on the APEX telescope 
(Giisten et al. 2006) at Llano Chajnantor, Chile. The maps are typically 40 arcseconds 
wide with a velocity channel width of 0.5 km s -1 . The analysis of these maps will be 



presented in a forthcoming paper, while here they will be only marginally used in Sect, 4.2 
for a comparison with the IRAC maps. 



3. Jet analysis 

3.1. Jet searching procedure 

YL09 have shown how the analysis of Spitzer data can be used to discover and 
characterize the emission from protostellar outflows. Such a characterization by means 
of the IRAC color-color plot relies on the presence of strong (0-0) and (1-0) H2 emission 
features in the IRAC bands (Neufeld & Yuan 2008), and makes it possible to define, on a 
[3.6] - [4.5] vs. [4.5] - [5.8] IRAC diagram, the region occupied by shocked H 2 as a function 
of the gas kinetic temperature and hydrogen density, in the range 1500 K < T < 4000 K 
and 10 2 cm -3 < n(H)< 10 5 cm -3 , respectively. 

In Figure [2] (adapted from Figure 1 of YL09), we show three regions (jet locus) 
identified by YL09 to contain jets: (1) the region corresponding to jets characterized by 
sufficiently high temperature and hydrogen density, that are unambiguously distinguished 
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Fig. 1. — Left panel - IRAC image of VMR-D at 4.5 /im (relative intensity is given in the 
bottom bar). White boxes represent the fields imaged with SofI at 2.12 /im (ISAAC for jet # 
15) where jets have been discovered (represented with red lines). Numbered labels refer to 
the jet identification numbers listed in Table[TJ Right panel - Map of the dust continuum at 
1.2 mm (adapted from Massi et al. 2007). The SofI fields (here shown in red) cover most of 
the brightest peaks. Intensity is given in the bottom bar in units of mJy/beam. 
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[4.5H5.8] 



Fig. 2.— Adapted from Figure 1 of YL09. IRAC [3.6] - [4.5] vs. [4.5] - [5.8] two colors 
plot indicating the region occupied by shocked H2 emission. Constant density (dotted) and 
temperature (solid) lines are indicated. Inside the diagram three regions are defined : zone 
1 - region encompassed by curves with 10 3 cm -3 < n(H) < 10 5 cm -3 and 2300 K < T < 
4000 K and having [4.5] - [5.8] < 0.8 mag; zone 2 - region possibly affected by contamination 
of Class 0/1 objects: i.e. within the curves with 10 2 cm~ 3 < n(R) < 10 5 cm -3 and 1500 
K < T < 2300 K and with [4.5] - [5.8] > 0.8; zone 3 - region empirically determined by 
Gutermuth et al. (2008) to contain outflows (dashed-dotted line, i.e. roughly the portion 
with [4.5] - [5.8] < 0.8 mag and [3.6] - [4.5] > 1.05 mag). The reddening vector, corresponding 
to Ay=30 mag, is depicted as well. 
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Fig. 3. — Left panel - Shock-mask around the location of jet # 4 of all the pixels having IRAC 
colors compatible with shocked H2. Colors (green, red and blue) are those corresponding 
to the three zones (1,2,3, respectively) depicted in Figure [2] Right panel - Same portion of 
the sky depicted in a three-color map obtained by combining the IRAC images in the three 
bands. The blue channel is used at 3.6 /im, the green one at 4.5 /mi, and the red one at 5.8 
/im. Note that, while the shock-mask evidences two possible jets in orthogonal directions, 
the three-color image confirms the presence only of the jet in the vertical direction (see text 
for more details). 
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from stellar objects; i.e. the region encompassed by curves with 10 3 cm 3 < n(H) < 10 5 
cm" 3 and 2300 K < T < 4000 K and having [4.5] - [5.8] < 0.8 mag, zone 1; (2) the region 
possibly affected by contamination of Class 0/1 objects: i.e. within the curves with 10 2 
cm" 3 < n(H) < 10 5 cm" 3 and 1500 K < T < 2300 K and with [4.5] - [5.8] > 0.8 mag, zone 
2; (3) the region empirically determined by Gutermuth et al. (2008) to contain outflow, i.e. 
roughly the portion with [4.5] - [5.8] < 0.8 mag and [3.6] - [4.5] > 1.05 mag, zone 3. 

Although the IRAC color-color plot represents a useful tool to search for jets, it has to 
be used with some caution, since it does not identify jets uniquely. Indeed, together with 
the contamination of zone 2 (see above), extinction effects can be significant in moving 
points inside and outside the jet locus. We show in Figure [2] the extinction vector for Ay 
= 30 mag; notably, even for extinction values lower than that, Class II sources may shift 
inside the jet locus. Moreover, possible emission of the CO v=l-0 fundamental transition at 
4.6 /im and [Fell] a 4 F-a 6 D lines between 4.4 /im and 4.9 /im can enhance the flux in band 2, 
moving the colors toward the top-left corner of the diagram. 

Therefore, as also suggested by YL09 themselves, the evaluation of the IRAC colors 
was complemented by the visual inspection of a three-color map of the region (constructed 
out of the 3.6 /im, 4.5 /mi, and 5.8 /im images) to examine the morphology of the emission 
in the pixels falling in the jet locus. 

In more detail, this part of the analysis was divided in two steps. First, the median 
subtracted IRAC images were converted into a color-map by applying the photometric 
corrections for extended sources suggested by the Spitzer Science Center website for low 
surface brightness measurements (correction factors of 0.91, 0.94 and 0.71 at 3.6, 4.5 and 
5.8 /im, respectively). As a result, we have obtained a mask for each of the three regions 
previously defined, which were then combined to obtain a false color representation (called 
'shock-mask') depicted for all the pixels having IRAC colors consistent with shocked H2. 
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As an example, we show in Figure[3j left panel, the enlargement of the shock-mask in the 
region of jet # 4, where the color code is: zone 1: green, zone 2: red, zone 3: blue. 

As a second step, we constructed a three-color image by combining the 3.6 /im (blue 
channel), the 4.5 /im (green channel) and the 5.8 /im image (red channel), with the aim 
to confirm (or not) by means of a visual inspection the possible shock-excited emission 
localized in the shock-mask. First, we discard all pixels that, although falling in the 
shock-mask, appear in the three-color image as isolated pixels or structures composed by 
less than 4-5 pixels which are likely due to artifacts or noise fluctuations. Then, we applied 
the following prescriptions to confirm groups of pixels in the shock-mask as real knots: 

- they must appear as extended structures. In SEC10, we define as point-like those 
sources with a 'shape parameter' Sh < 0.7, where Sh=a 2 (observed) - <r 2 (PSF), is 
the difference between the squares of the half-width of the best-fitting function and 
the point-spread function. Here, we request that knots must have Sh > 0.7. This 
prescription, although implying that possible point-like knots are discarded, has 
however the advantage of easily individuate protostellar sources in the jet locus of the 
color-color diagram. 

- preferably, knots should belong to a chain (which can have a straight-, curved-, or 
S-shape). They must be detected at signal-to-noise level greater than 3 in at least 
one band. 

- isolated knots are considered only if detected at signal-to-noise level greater than 3 in 
at least two bands. 

In Figure[3j right panel, we show the enlargement of the three-color image in the region 
of jet # 4. This is an illustrative case of the importance of the visual inspection described 
above: while the shock-mask in the left panel evidences two possible jets in orthogonal 
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directions, the visual inspection of the three-color image confirms the presence only of the 
jet in the vertical direction, since the pixels in the horizontal direction correspond to image 
artifacts, to point-like sources, or to pixels whose photometry is strongly contaminated by 
the emission of the central cluster. 

3.2. Jet description 

With the procedure described in the previous section, we have found a total of 15 jets 
inside the area observed with IRAC, as well as some sparse knots that do not show any clear 
alignment. Noticeably, we re-confirm with IRAC all the jets that were already detected in 
the near-infrared, except jets # 7 and # 14 that represent new discoveries. This is depicted 
in Figure[TJ where it is clear that almost all the Spitzer jets (red lines) are located within 
a field observed in the H 2 2.12 /iin line. Indeed, such result is not unexpected, since H 2 
Spitzer and 2.12 /im lines, coming from levels with similar excitation energy, probe similar 
excitation conditions. Conversely, some knots (within the jets) detected in the 1-0 S(l) line 
with SofI do not appear in any Spitzer image. In the case of jet # 15, this occurs because 
the most internal knots (e.g. knots D and F reported by Giannini et al. 2005) are confused 
in the Spitzer images in the strong diffuse emission of the protostellar cluster where the 
exciting source is located. In other cases of faint knots (at 2.12 /im), the lack of a Spitzer 
counterpart can be ascribed to the poorer sensitivity at longer wavelengths. For example, 
the 3 a sensitivity limit at 2.12 //m of ~ 0.01 mJy is around a factor of 3 better than the 
sensitivities at 3.6 /im and 4.5 /im. 

The very nice coincidence between IRAC and SofI jet detections leads to three main 
conclusions: i) the approach adopted to discover jets in the IRAC images can be considered 
as validated, ii) H2 2.12 /im surveys allow to catch a relevant fraction of jets and hence 
longer wavelengths observations are not mandatory to overcome extinction problems, and, 
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Fig. 4. — Jet 1: Left panel - Three colors (3.6 /im : blue - 4.5 /im : green - 5.8 /im : red) 
IRAC image. H2 knots locations are highlighted with boxes and labeled as in Table[2j The 
candidate exciting source(s) is indicated with a circle. Right panel - Three colors (3.6 /im 
: blue - 4.5 /xm : green - 24 /im : red) IRAC/MIPS image where the contours of the dust 
emission at 1.2 mm (Massi et al. 2007) are superposed in steps of 4 a starting from a level 
of 50 mJy/beam. The candidate exciting source(s) is (are) indicated with a circle. 
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Fig. 5. — As in Figure [4] for jet 2. 
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Fig. 6. — As in Figure [4] for jet 3. 
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Fig. 7. — As in Figure [4] for jet 4. 
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Fig. 8. — As in Figure [4] for jet 5. 
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Fig. 9. — As in Figure [4] for jet 6. 
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Fig. 10. — As in Figure [4] for jet 7. 
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Fig. 11. — As in Figure [4] for jet 8. 
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Fig. 12. — As in Figure [4] for jet 9. 
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Fig. 13. — As in Figure [4] for jet 10. 
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Fig. 14. — As in Figure [4] for jet 11. 
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Fig. 15. — As in Figure [| for jet 12. 
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Fig. 16. — As in Figure [4] for jet 13. 




Fig. 17. — As in Figure [4]for jet 14. 
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Fig. 18. — As in Figure [4] for jet 15. 
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%%%) the absence of jets in the parts of the IRAC map not covered also in the near-infrared 
strongly indicates that jets can be preferentially found in proximity of the dust peaks, 
namely where the near- infrared images were originally taken (see also De Luca et al. 2007). 



The 15 Spitzer jets are depicted in Figures|4] to 18 and are described in Appendix |Aj 



In each of these figures, we show in the left panel a composite IRAC 3-color image and in 
the right panel a composite IRAC-MIPS image with overlaid the contours at 1.2 mm (these 
latter in steps of 4 a starting from a 4 a level of 50 mJy/beam). In these figures, jet knots 
are marked by small boxes while the candidate exciting source, determined as described in 
Sect.[4j is marked by a small circle. 

A summary of the jet properties is given in Table [TJ where for each jet we give in 
columns 1-4 the identification numbeiQ the number of the detected knots, the projected 
length, and the position angle. As anticipated, all the jets were originally discovered in the 
H 2 2.12 /im line, the only exceptions being represented by jet #7, # 14 and by other few 
knots (K1-K4), which were not targeted in the near-infrared. In the last column, we give 
the references where these jets have been reported. The projected jet lengths and position 
angles are in agreement with those given in the literature, the only exception being the 
length of jet # 6, which was originally supposed to have a double length with respect to 
that reported in Table[l] (Giannini et al. 2007). This hypothesis was done since the jet 
appears one-sided in the 2.12 /im image, and the lack of a counter-jet was ascribed to the 
limited field of view of Sofl. Based on the IRAC image (see Fig.[9]) we now confirm the lack 
of any additional knot symmetric to the candidate ES. As a general note, jets # 2 ,3 ,7, and 
11 are rather compact (length lower than 0.15 pc), jets # 4 ,14, and 15 extend for more 
than 0.4-0.5 pc, while all the others have intermediate lengths. Length of jets #5 and # 13 
is very uncertain, since it is not clear whether or not the discovered knots can be effectively 



x The jets identification number follows our internal classification. 
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attributable to the same jet (see also the discussion in AppendixjAJ) . Regarding the position 
angle (PA, computed from the North to the East), notably around half of the jets are 
oriented 160° ± 15° : we will comment this feature in Sect. |5.3l 



3.3. Jet photometry 

To obtain the photometry of the detected knots we have applied the task polyphot 
within the Image Reduction and Analysis Facility (IRAF) package]^] which is suited for 
evaluating the emission of extended and irregular sources. For each knot (both in the IRAC 
bands and in the narrow-band at 2.12 /xm), we have first defined a polygon with sides 
tangent to the 3 a contour of the emission and then summed up the flux of all the pixels 
inside it. The sky background was estimated locally and then subtracted from the flux 
computed inside the polygon. As shown by Neufeld & Yuan (2008), the photometry in the 
IRAC bands are dominated by bright H 2 lines, and in particular by the v—1-0 0(5)-0(7) 
in band 1 (3.6/xm), t>=0-0 S(9) in band 2 (4.5 /xm), v—0-0 S(7)-S(6) in band 3 (5.8/jm), 
and v=0-0 S(5)-S(4) in band 4 (8.0/jm). Minor contributions are expected from PAH, CO 
and atomic lines which, however, cannot be easily estimated (Takami et al. 2010) and thus 
are not taken into account by us in the photometry. Similarly, we ascribe the near-infrared 
narrow-band filter photometry totally to H2 1-0 S(l) line emission, since no continuum 
emission in the fT-band is usually seen in the spectra of protostellar jets (e.g. Giannini et al. 
2004, Nisini et al. 2002). Table[2] lists the coordinates and the size of each knot, evaluated 
on the 4.5 fim image, along with the photometry between 2.12 /im and 8.0 fim. The last two 
columns list the IRAC luminosity (£irac> obtained by summing up the luminosities of all 
the IRAC bands) and 2.12/xm luminosity (Li_os(i))j both computed assuming a distance to 



2 Available at 



http:/ /iraf.noao.edu/ 
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VMR-D of 700 pc (Liseau et al. 1992). 

4. Exciting sources analysis 

4.1. Exciting sources identification 

As a second goal of our study, we have used the IRAC and MIPS images to search for 
the ES using the following procedure. 

- On the IRAC image at 4.5 /im we have defined a rectangle oriented as the jet itself 
and having the long side equal to n x the jet length, where the value of n, in the 
range 1-3, is selected for each jet to not exceed the parsec scale. The short side, 
orthogonal to the jet direction, is defined case-by-case so that all the knots lie inside 
the rectangle. 

- We then searched the IRAC catalog (SEC10) for all the point-like sources inside the 
rectangle, selecting those whose position is compatible with the jet shape (straight- 
curved- or S-shape). 

- Sources are then selected if they have a spectral index (computed by using the 
available measurements between K s and 8 /im) compatible (within the upper limits) 
with a YSO of Class II, I or flat. To Class I or flat sources we give higher priority than 
to nearby Class II sources. Analogously, we have searched the MIPS catalog to select 
possible sources not present in the IRAC catalog because not detected at wavelengths 
shorter than 24 /zm. 



The results of this search are presented in Table[3j We remark that, given the crowding 
of YSOs in the investigated fields, we can not consider our search exhaustive, especially for 
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the presence of many Class II sources in the defined rectangles. Consequently, the sources 
listed in Table[3]have to be considered as ES candidates and hence as interesting objects to 
be further studied. Noticeably, however, we find only one best candidate in all cases, apart 
from jets # 5 and # 14, for which two possible candidates have been selected (labeled with 
a and b according with a priority order, see AppendixlAl) . 



Associations with counterparts at different wavelengths (from 1.2 /im to 1.2 mm) 
are reported, with the corresponding identification number assigned in each catalog. 
Associations between IRAC, MIPS and SIMBA sources have been already done by SEC10 
and reported in their photometric catalog; here we also consider possible associations with 
2MASS, WISE and BLAST sources. These have been matched to the IRAC source list by 
using a circle whose radius is the sum of the positional uncertainties associated with each 
pair of instruments (reported in the caption of Table[3]). The ES of jets #2, #3, # 11, # 13 
are very barely or not detected at wavelengths shorter than 24 /im; this confirms the results 
presented in Giannini et al. 2007 (for jets #2, #3) who hypothesized them as very young 
YSOs. 

All except one sources have a MIPS 24 /mi counterpart and 8 (10) are associated with 
a 250 (1.2 mm) peak, while only 5 objects are seen in the 2MASS bands. The source 
exciting jet # 15 has been already the subject of a dedicated study (Giannini et al. 2005). 
It is a member of an infrared cluster that appears completely unresolved and saturated in 
the Spitzer bands. For this reason that object (also known as IRS 17-# 40-3) will not be 
commented further on in the present paper. 



In Table|4]we list the fluxes (in mJy) of the ES, which will be used in Sect. 5.1 to derive 
their SED. As far as the WISE data are concerned, we give only the fluxes at 12 /im (band 
3). Those in bands 1,2,4 (at 3.4 jum, 4.6 /im, and 22 jum) do not differ from the Spitzer 
ones (at 3.6 /im, 4.5 /an, and 24 /im) by more than 30%, which does not imply significant 
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modifications in the SED global shape. Two sources (# 4 and # 14a) are saturated at 
24 /iin, hence the WISE 4 band flux has been used for their SED fitting. Finally, in the last 
column of Table[4] we list the source spectral index, computed for sources with at least four 
valid fluxes between 2.16 /im and 12 /iin. 



4.2. Evaluation of the jet extinction 

As a further step toward identifying the driving sources, we have investigated whether 
the selected candidate ES of each jet is located at the center of two distinct lobes traveling in 
opposite directions, as expected for bipolar jets. Since we lack high-resolution spectroscopy 
observations, which would have allowed us to directly measure the mutual shift of the 
emission lines, we have adopted an indirect way to identify (if any) the blue- and red-shifted 
lobes. This method is based on the comparison of the IRAC total luminosity (Xirac) 
with the luminosity of the 1-0 S(l) line (L 1 _ 5(i)) ) both listed in the last two columns of 
Table[2j Such comparison helps in identifying the blue- and red-shifted lobes under some 
conditions: i) the visual extinction increases in the direction of the red lobe (as expected if 
this coincides with the receding part of the jet and if this latter is not in the sky-plane); ii) 
Li_os(i)/£irac weakly depends on the temperature and density in a given knot. 



To verify the condition ii), we have estimated the ratio -^i-os(i)/-^irac f° r temperatures 
between 1000 K and 3000 K, in the Local Thermal Equilibrium (LTE) approximation. As 



shown in Figure 19, such ratio remains fairly constant in the range ~ 0.2-0.4 (note that this 
is generally not true if the ratio of Li_ Q s(i) with the luminosity in a single IRAC band is 
taken). Hence, a difference of -^i-os(i)/-^irac of more than a factor of ~ 2 in symmetrically 
located knots (with respect to the putative exciting source) can be primarily ascribed to a 
difference in the local extinction (A Ay), since this latter affects the near-infrared and the 
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mid-infrared fluxes in a differential way. (We note, however, that the assumption ii) is too 
coarse to allow the absolute determination of Ay). 

In Table[5]we list for each jet the pairs of symmetric (groups of) knots for which we 
have evaluated A Ay, computed by assuming the Rieke & Lebofsky (1985) extinction law. 
To minimize the uncertainties, this has been obtained from the average £i_os(i)/£irac i n 
each side of the jet. Adding the photometric uncertainties to the allowed variations in 
L 1 _ 5( 1 )/L IR Ac due to the physical conditions in the knots, we consider as significant only 
AAy larger than 2 mag, which roughly correspond to a difference of Li_os(i) /£ibac °f 
more than 4. Summarizing the results, our objects are divided in three groups, in which we 
measure a AAy with a decreasing level of significance. In particular we get: AAy >5 mag 
in five jets (# 1, 2, 4, 8, 11), 2 mag < AAy < 5 mag in two jets (# 3, 12) and AAy < 2 mag 
in four jets 5, 9, 10, 15). For jet #13 we measure AAy both considering knots 3 and 
4 (AAy= 10 mag) and the groups of knots 1-2 and 5-6 (AAy= 4 mag). In both cases all 
the knots south (north) the ES are red- and blue-shifted, respectively. This result slightly 
favors the hypothesis that all the knots are part of the same jet (see also Appendix[Aj . For 
the remaining jets we are not able to apply the method either because we can not identify 
knots symmetrically located (jet # 6) or because measurements of the 2.12 /im line do not 
exist (jets # 7, 14). In particular, for jet # 14 (as well as for jet # 5, where AAy is not 
significant) we can not use the 'extinction method' to identify the exciting source between 
the two candidates labeled with a and b in TableEl 

Remarkably, the few direct velocity measurements existing in the literature for the 
VMR-D jets are all in agreement with the results obtained with the above 'extinction 
method'. These are the high-resolution spectroscopy observations of the 2.12 /im line of 
jet# 15 (Giannini et al. 2005), the 12 CO(1-0), 13 CO(2-l) maps in the vicinity of the IRAS 
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sources and overlapping jets # 4, 8, and 15 (Elia et al. 2007), and the small 12 CO(3-2) 
maps obtained with APEX for jets # 1, 2, and 3 (see as an example the map of jet # 3 



shown in Figure 20). Another 12 CO(3-2) map (although incomplete) of jet # 5 should favor 
source 5b as the exciting source, since all the parts of the outflow overlapping knots 2-7 
appear blue-shifted. 



5. Discussion 
5.1. SEDs fits 

Once the observed fluxes have been determined to construct the SEDs for our candidate 
ES, we used the radiative transfer model developed by Whitney et al. 2003a, Whitney et 
al. 2003b to obtain insight on the physical properties of these objects. Actually, the use of 
this model has been greatly facilitated by Robitaille et al. 2006, Robitaille et al. 2007, who 
provided a grid of 200,000 computed models and an efficient online fitting procedurd^to 
help the community in the SED fitting work. In the best-fit procedure we set as constraints 
the extinction between mag and 200 mag and the distance to the VMR-D within the 
range 500 pc - 1000 pc (Liseau et al. 1992). We fit the SED of sources with at least six 
photometric detections in Table|4j In addition, a single grey-body function was used to fit 
the SED of the sources not observed at wavelengths shorter than 24 fim, but with at least 3 
data points (# 2, # 3, # 13). We used different-sized apertures with radius corresponding 
to 3 a of the instrument beam and verified that even adopting larger apertures the same 



models are selected. The results of the best-fits 'a la Robitaille' are depicted in Figure[21 
Here we show as upper limits also the BLAST (350-500 /mi) data points, although not 
considered in the fitting procedure. This choice is motivated by the fact that multiple 



3 Available at 



http:/ /caravan. astro, wisc.edu/protostars/sedfitter.php 
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sources lie in the large BLAST beam. The main output parameters are reported in Table[6j 
where we give both the values of the best-fit model (in boldface) and those of the model 
whose x 2 is 50% higher than the minimum value. These parameters are the distance to 
the source (col. 2), the disk inclination (between 0° and 90° with respect to the sky-plane, 
col.3), the interstellar and circumstellar extinction (cols. 4, 5), the stellar temperature and 
mass (cols. 6, 7), the outer envelope radius (col. 8) along with its mass (col. 9), and the total 
(stellar plus accretion) luminosity (col. 10). The general picture emerging from the SEDs 
modeling is that most of the ES are low-mass and relatively evolved protostars. The stellar 
temperatures (2900 K < T* < 4400 K) and masses (typically M ps 0.05-1 M ) are those 
typical of embedded (Ay ~ 10 - 60 mag) T Tauri stars of spectral type between K0 and 
M0. The reservoir of material of the surrounding envelope (M env typically between 10 -2 - a 
few M ) is not enough to enhance significantly the final stellar mass. Possible exceptions 
to this general picture are : i) sources # 1 and # 4, whose mass (already piled up onto 
the star or still in the envelope) appears potentially sufficient to form an intermediate-mass 
object. Given the crowding of the regions where these sources are located (close to infrared 
clusters, see also Table[T]), we can not exclude however that they are binaries or multiple 
objects unresolved at the IRAC angular resolution; ii) sources # 7 and # 14a, for which 
the fitted M env and R env are very uncertain, since their SED is not measured beyond 
A=70 /im. Finally, we note that the disk inclinations of sources # 1, # 4, # 8, for which 
the correspondent jets have two mutually shifted lobes (i.e. for which AAy is > 5 mag, 



see Sect, 4. 2), are 20°-40°, 55°-65°, 20°. In the hypothesis that jet and disk are orthogonal, 
these inclinations confirm that these jets lie outside the sky-plane. 
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5.2. Jets vs. exciting sources properties 



In Figure 22 the projected jet length (in pc) is displayed as a function of the far-IR 
color [24]-[70jjof the candidate ES. This plot is based only on sources selected as first 
priority ES (see also Appendix[A| . The idea underlying such plot is to investigate whether 
or not a relationship exists between the two considered parameters within an evolutionary 
framework. Indeed, the youngest jets are expected to be the shortest, in the assumptions 
that all the jets travel at comparable velocities and that their inclinations with respect to 
the sky-plane are not dramatically different. In this sense the following analysis has to be 
considered in a statistical sense, and no conclusions can be drawn for the individual sources. 
The [24]- [70] color is typically used clS cLll cl ge indicator, because the youngest sources 



are also the coldest and reddest ones. Figure 22 shows an evident correlation (regression 



coefficient > 0.8) which well agrees with the expectations. Noticeably, the same kind of plot 
constructed with IRAC (instead of MIPS) colors or with the spectral index between 2 fim 
and 12 /zm does not show any correlation: bluer colors are, in fact, more sensitive both to 
different circumstellar morphologies (e.g. disk inclination) and to extinction effects. On the 
other hand, the [24]- [70] color is more sensitive to the circumstellar envelope temperature 
and almost unaffected by extinction, since it changes of only 0.2 mag for Ay 100 mag. 
Hence, the observed correlation is compatible with an evolutionary interpretation, and can 
not be related to a mere inclination effect: if this latter were responsible for the shortening 
of the jet projected length, the bluest (warmest) portions of the circumstellar disk material 
would be exposed to the observer, hence an anti-correlation would emerge between the jet 
projected length and the [24]- [70] color. 



In Figure 23 we plot the IRAC line cooling as a function of the ES bolometric luminosity 



4 The MIPS [24]-[70] color is defined as 2.51og(F 70 /F 24 x F0 24 /F0 70 ) where F 70 ,F 24 are 
the fluxes at 24, 70 fim, respectively and -FO70, -FO24 the zero magnitude fluxes. 
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(considering only ES of higher priority). Our aim is to investigate the balance between the 
energy irradiated by the driving protostar and that emitted by the associated jet. In that 
Figure, the horizontal dotted lines are the range of bolometric luminosity given in Table|6j 
with the value associated with the best-fit model indicated with a dot. The reported values 
of Lirac are those given in Table[2] without any extinction correction; for reference, these 
values increase by a factor of 0.3 and 1.3 (in logaritmic scale) for Ay = 10 mag and 30 mag, 
respectively. Instead, the uncertainties on Lirac associated with the photometric errors are 
not significant. A rough linear correlation (in the log- log scale) is found between the two 
quantities (regression coefficient of 0.68), with an intercept value located at Lirac ~ 10~ L96 
L ). This value can be easily converted in the total H2 cooling (Lh 2 ) by estimating the 
ratio Lh 2 /£irac under LTE conditions: we find that this ratio is ~ 4 at temperatures of 
2000-3000 Kelvin. Hence, the above relation translates into Lh 2 ~ 10~ 156 L . This value 
is in good agreement with the result (Lh 2 ~ 10~ 1,4 L ) found by Caratti o Garatti et al. 
(2006) in a sample of Class 0/1 sources. In that work no significant differences in terms 
of efficiency were found between the two classes, likely because of the large uncertainties 
associated with Lh 2 , which was largely derived on the base of near-infrared observations. 
These differences are indeed expected in the commonly accepted scenario in which a tight 
relationship exists between mass-loss and mass accretion rate, with a declining of the latter 



with the source evolution (e.g. Bontemps et al. 1996), In Figure 23 a certain degree of 
segregation between the points above (sources # 1, 3, 4, 6, 9, and 13) and below the 
straight line (sources # 2, 7, 8, 10, and 14a) is indeed recognizable, the formal fits through 
the two groups having regression coefficients of 0.97 and 0.99, respectively. We estimate 
that the two groups can mix each other only if their differential extinction is (on average) 
larger than 15-20 mag. Even a more clear segregation is found by using far-infrared lines. 
For example, the ratio between Lfir (i.e. the cooling rate due to the far-infrared lines) with 
-^boi gives values spanning from 10~ 2 to 10~ 4 going from Class to Class II sources (Nisini 
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et al. 2002, Lorenzetti et al. 2002a). 

5.3. Large-Scale Properties 

The circumstellar accretion disks are expected to have their symmetry axis parallel to 
the local magnetic field (Shu et al. 1995). At larger scales, the magnetic field is expected to 
have a role only in regions located in the Galactic plane and far from the center, i.e. where 
it is oriented parallel to the Galactic disk (150°-160°, Stephens et al. 2011). Indeed, jets 
with orientations parallel to the large-scale magnetic field have been already found both in 
the Galactic plane (Hodapp 1984, 1987; Scarrott et al. 1992; Lorenzetti et al. 2002b) and 
outside of it (in Taurus - Menard & Duchene 2004). 

Given its location, VMR-D represents a well suited case to test the influence of the 
Galactic magnetic field on the jet orientation. Remarkably, out of 15 jets, 8 (# 1, 4, 6, 8, 
9, 11, 13, 15) have a position angle of ~ 160°± 15° (see Table[T]). For completeness, also 
the jet from IRS8 (located in VMR, as well) presents an orientation of ~ 153° (Lorenzetti 
et al. 2002b). Hence, the present survey suggests that the Galactic magnetic field may be 
important during the early phases of the pre-main sequence evolution. However, a selection 
effect could bias our sample and make this inference weaker: since the influence of the 
Galactic magnetic field tends to decrease toward the inside of the cloud, we expect that jets 
aligned with it are located preferentially in the external layers. Therefore, these jets are 
likely the least embedded and therefore those most easily detected. 

5.4. Comparison with other H2-jet surveys 

The present H 2 -jet survey covers a rather small region (about 1.2 deg 2 ), but it has 
the advantage of being substantially unbiased and therefore complete at the Spitzer 
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sensitivity. Hence it is worthwhile to compare our results with those stemming from other 
similar surveys that, having both investigated larger areas and detected a larger number of 
sources, offer a more robust statistical significance in tracing the early phases of the stellar 
formation. Indeed, the peculiarity of our survey is to study a region inside a GMC, which 
is the typical site of the star formation, but rather different from other nearer dark clouds 
where star formation occurs in a usually less turbulent environment. 

There are currently two main projects regarding 2.12 /zm H 2 jet emission: (i) the survey 
UWISH2 (UKIRT Widefield Infrared Survey of H 2 , Froebrich et al. 2011); and (ii) MHOs 
(Molecular Hydrogen Emission-Line Object^]), namely the general catalog of H 2 infrared 
(1-2.2 fim) emission from young outflowing objects (Davis et al. 2010). On both papers 
references are given to the numerous 2.12 /im H 2 imaging surveys of different regions carried 
out in the past years. Since UWISH2 will eventually provide an unbiased and complete 
census of the 1-0S(1) ro-vibrational H 2 line emission along the Galactic Plane, beside star 
forming and HII regions, it will include also evolved (AGB) stars, supernova remnants 
and fluorescent H 2 emission. Conversely, MHOs catalog is a fundamental database only of 
YSOs jets and HH objects that present molecular emission. Some of the jets described in 
the present paper have been already mentioned in a very preliminary form (Giannini et 
al. 2007), and, as such, are included in MHOs. Finally, an example rather similar to our 
VMR-D survey (as far as its GMC nature, size of the investigated area, distance from the 
Sun are concerned) is represented by the extensive Orion A survey (Stanke et al. 2002). 
Some results from the literature are the following: 

UWISH2 (a) A potential source candidate can be assigned to about half of the jets; (b) 

typically, the flows are clustered in group of 3-5 objects, within a radius of 5 pc; (c) 
these groups are separated by about 5 pc; (d) 2/3 of the surveyed area is devoid of 



5 Available at 



http://www.jach.hawaii.edu/UKIRT/MHCat/ 
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jets; (e) a large range of apparent lengths is found (from 4 to 130 arcsec); (f) a small 
fraction (about 10%) are "parsec scale" jet; (g) there is a marginal (less than 3a) 
over-abundance of flow position angles roughly perpendicular to the Galactic Plane. 
All of these are preliminary considerations based on the Serpens and Aquila survey 
(Ioannidis & Froebrich, 2012). 

MHOs (a) in some regions multiple knots and bow shocks radiate in many directions from 
a tight cluster of young stars and the relationship between these objects is often 
unclear, mainly in massive star forming clusters; (b) jets are found in both low- and 
high-mass star forming regions. 

Orion A (a) the first outcome of this survey is the large number of detected jets distributed in a 
more isolated way with respect to Orion B that presents a more clustered distribution; 
(b) a large variety of morphologies have been detected, ranging from large, extended, 
filamentary features to compact or even unresolved knots. This circumstance makes 
the identification of flows not always an objective process; (c) for a reasonable 
fraction of the jets a candidate exciting source is either obvious or suggested; (d) few 
jets show some degree of symmetry; (e) very few jets are well collimated, virtually 
unresolved orthogonally to the flow direction; (f) apparently, there is not a preferred 
flow orientation. 

The results of our VMR-D survey differ in several respects from the previous surveys' 
findings, mainly because these are based only on near-IR data. In the majority of cases 
(13 out of 15) a convincing exciting source has been found: this circumstance is related to 
the increased capability offered by the longer wavelengths coverage to identify young and 
embedded objects. Practically, all the VMR-D jets display a high degree of collimation 
and most of them are composed of knots arranged in symmetrical position. Preserving 
the collimation and a certain symmetry over a large part of the GMC is a key-feature 
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that makes easier the identification of a candidate driving source and it could signal the 
presence of a relatively quiet environment, not too much affected by (local) turbulence. In 
this sense, VMR-D could be different from both the more active VMR-A,-B,-C clouds and 
Orion GMC; indeed, this picture is also supported by independent evidences as the lack in 
VMR-D of O-B associations (Liseau et al. 1992) and the overall gas dynamics (Olmi et al. 
2010). 

6. Summary 

This paper presents a systematic search of molecular jets from young protostars in 
the star formation region Vela, cloud D. Our color-color analysis in the four bands of 
Spitzer-YRAC has led to the following results: 

- 15 jets have been detected in an area of about 1.2 square degrees, as well as a few 
sparse knots not displaying a clear alignment. All except two jets were previously 
known since already imaged in the H 2 2.12 /zm line. The absence of jets in the part of 
the IRAC map not covered also in the near-infrared strongly indicates that jets can 
be preferentially found in proximity of the dust peaks, since these are the locations 
where the 2.12 /xm images have been taken. 

- Using the available catalog of YSOs in VMR-D we searched the jet exciting sources. 
For each jet, we selected all nearby Class I, flat, and Class II sources, and selected 
those at the jet center and aligned with it. In this way we are able to find a 
best-candidate exciting source in all except two jets, for which we remain with two 
possible alternatives. To validate the jet-ES association, we have investigated whether 
a differential extinction exists between the two (hypothetical) lobes. Through an 
analysis of L H2 /£irac; we find a significant extinction gradient of A Ay > 5 mag in 
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five cases. 

- A search for the ES counterparts at wavelengths between 1.2 /xm and 1.2 mm has 
been done. All except one source have a 24 /xm counterpart, and 10 are associated 
with a 1.2 mm peak. Four sources are very barely or not observed at wavelengths 
shorter than 24 /xm, suggesting they are very young protostars; three of them are also 
associated with the most compact jets (I < 0.1 pc). 

- The SEDs of the ES have been constructed based on all available photometric data 
between 1.2 /xm and 1.2 mm. We used the Robitaille et al. (2006, 2007) model to 
derive the main source parameters, such as stellar temperature and mass, envelope 
mass and radius, and total bolometric luminosity. On average, the ES are low-mass 
and relatively evolved protostars, with the possible exception of few objects located 
close to infrared clusters, which could be intermediate-mass protostars as well as 
binaries or multiple objects not resolved at the Spitzer angular resolution. 

- A significant correlation is found between the projected jet length and the [24] - [70] 
color, which is consistent with an evolutionary scenario according to which shorter 
jets are associated with younger sources. We remark, however, that this trend has to 
be considered just in a statistical sense, and that no conclusions can be drawn for the 
individual sources, given the strong assumptions (e.g. jet inclination) under which 
the correlation is derived. The same trend is not recognizable if the diagram is drawn 
with shorter wavelengths colors, whose values are significantly affected by geometrical 
and extinction effects. 

- A coarse correlation is found between IRAC line cooling and the ES bolometric 
luminosity, which is in agreement with previous literature determinations. It is 
suggested that the sources segregate in two distinct groups, to which pertain a 
different efficiency in the line cooling. This trend is in agreement with the the scenario 
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in which mass loss and mass accretion are tightly related phenomena, with mass 
accretion rates decreasing with time. 

- About half of the jets present a position angle which may indicate a certain influence 
by the Galactic magnetic field, at variance with other star forming regions not located 
exactly in the Galactic plane. However, a selection bias could affect our sample. 
Indeed, since the influence of the Galactic magnetic field tends to decrease toward the 
inside of the cloud, we expect that jets aligned with it are located preferentially in the 
external layers, where they are most easily detected. 

- Unlike other similar studies conducted on other star forming regions, the association 
between jet and ES has been more successful in VMR-D. Together with the increased 
capabilities offered by the longer wavelengths at which VMR-D is observed, this 
circumstance is likely favored by the high degree of collimation displayed by the jets. 
This, in turn, reflects the presence of a relatively quiet environment, for example in 
comparison with the other Vela clouds where star formation is going on. 



7. Acknowledgements 

This paper is dedicated to the memory of Alberto Salama who prematurely passed 
away during the development of the work. 

It is based on observations made with the Spitzer Space Telescope, which is operated 
by the Jet Propulsion Laboratory, California Institute of Technology under a contract with 
NASA. 



-40 - 
REFERENCES 

Black, J.H., & van Dishoeck, E.F. 1987, ApJ, 322, 412 

Bontemps. S., Andre, P., Terebey, S., & Cabrit, S. 1996, A&A, 311, 858 

Caratti o Garatti, A., Giannini, T., Nisini, B., & Lorenzetti, D. 2006, A&A, 449, 1077 

Caswell, J. L., & Haynes, R. F. 1987, A&A, 171, 261 

Cuby, J. C, Lidman, C, Johnson, R., Jaunsen, A., Moutou, C. 2004, ISAAC User Manual 

Cyganowski, C. J., Brogan, C. L., Hunter, T. R., Churchwell, E., Zhang, Q. 2011, ApJ, 729, 
124 

Davis, C.J., Froebrich, D., Stanke, T., et al. 2009, A&A, 496, 153 

Davis, C.J., Gell, R., Khanzadyan, T., Smith, M.D., & Jenness, T. 2010, A&A, 511, A24 

Davis, C.J., Scholtz, P., D., Lucas, P., Smith, M.D. & Adamson, A. 2008, MNRAS, 387, 954 

De Luca, M., Giannini, T., Lorenzetti, D., et al. 2007, A&A, 474, 863 

Elia, D., Massi, F., Strafella, F., et al. 2007, ApJ, 655, 316 

Fazio, G.G., Hora, J.L., Allen, L.E. et al. 2004, ApJS, 154, 10 

Froebrich, D., Davis, C. J., Ioannidis, G., et al. 2011, MNRAS, 413, 480 

Giannini, T., McCoey, C, Caratti o Garatti, A., et al. 2004, A&A, 419, 999 

Giannini, T., Lorenzetti, D., De Luca, M., et al. 2007, ApJ, 671, 470 

Giannini, T., Massi, F., Podio, L., et al. 2005, A&A, 433, 941 

Glisten, R., Nyman, L. A., Schilke, P., et al. 2006, A&A, 454, 13 



-41 - 

Gutermuth, R.A., Myers, P.C., Megeath, S.T., et al. 2008, ApJ, 674, 336 

Hodapp, K.-W. 1984, A&A, 141, 255 

Hodapp, K.-W. 1987, ApJ, 319, 842 

Hollenbach, D., & McKee, C.F. 1989, ApJ, 342, 306 

Kaufman, M.J., & Neufeld, D.A. 1996, ApJ, 456, 611 

Ioannidis, G., & Froebrich, D. 2012, MNRAS, 421, 3257 

Lee, H.-T., Takami, M., Duan H.-Y., et al. 2012, ApJS, 200, 2 

Lidmann, C, Cuby, J.-G:, Vanzi, L., et al. 2006, SofI User's Manual 

Liseau, R., Lorenzetti, D., Nisini, B. Spinoglio, L. & Moneti, A. 1992 A&A, 265, 577 

Lorenzetti, D., Giannini, T., Nisini, B., et al. 2002a, A&A, 395, 637 

Lorenzetti, D., Giannini, T., Vitali, F., Massi, F. & Nisini B. 2002b, ApJ, 564, 839 

Massi, F., De Luca, M., Elia, D., et al. 2007, A&A, 466, 1013 

Massi, F., Giannini, T., Lorenzetti, D., et al. 1999, A&AS, 136, 471 

Massi, F., Lorenzetti, D., & Giannini, T. 2003, A&A, 399, 147 

Massi, F., Lorenzetti, D., Giannini, T., & Vitali, F. 2000, A&A, 353, 598 

Menard, F., & Duchene, G. 2004, Ap&SS, 292, 419 

Murphy, D.C. & May, J. 1991, A&A, 247, 202 

Netterfield, C. B. Ade, P. A. R., Bock, J. J., et al. 2009, ApJ, 707, 1824 
Neufeld, D.A., & Yuan, Y. 2008, ApJ, 678, 974 



Nisini, B., Caratti o Garatti, A., Giannini, T., et al. 2002, A&A, 393, 1035 
Nisini, B., Giannini, T., & Lorenzetti, D. 2002, ApJ, 574, 246 
Noriega-Crespo, A., Morris, P., Marleau, F. R., et al. 2004, ApJS, 154, 352 
Ohlendorf, H., Preibisch, T., Gaczkowski, B., et al. 2012, A&A, 540, 81 
Olmi, L., Angles-Alcazar, D., De Luca, M., et al. 2010, ApJ, 723, 1065 
Pascale, E., Ade, P.A.R., Bock, J.J., et al. 2008, ApJ, 1081, 400 
Rieke, G. H., & Lebofsky, M. J. 1985, ApJ, 288, 618 

Robitaille, T.P, Withney, B., Indebetouw, R., Wood, K., & Denzmore, P. 2006, ApJS, 167, 
256 

Robitaille, T.P, Withney, B., Indebetouw, R., & Wood, K. 2007, ApJS, 169, 328 

Scarrott, S.M., Draper, P.W., Rolph, CD., & Wolstencroft, R.D. 1992, MNRAS, 257, 485 

Shu, F.H., Najita, J., Ostriker, E.C., & Shang, H. 1995, ApJ, 455, L155 

Smith, M.D. 1995, A&A, 296, 789S 

Smith, M.D. & Rosen, A. 2005, MNRAS, 357, 1370 

Smith, H.A., Hora, J.L., Marengo, M., & Pipher, J.L. 2006, ApJ, 645, 1264 

Stanke, T., McCaughrean, M.J., & Zinnecker, H.D. 1998, A&A, 332, 307 

Stanke, T., McCaughrean, M.J., & Zinnecker, H. 2002, A&A, 392, 23 

Stephens, I.W., Looney, L.W., Dowell, CD., Vaillancourt, J.E., & Tassis, K. 2011, ApJ, 
728, 99 



Strafella, F., Elia, D., Campeggio, L., et al. 2010, ApJ, 719, 9 (SEC10) 

Skrutskie, M.F., Cutri, R.M., Stiening, R., et al. 2006, AJ, 131, 1163 

Takami, M., Karr, J.L., Nisini, B., et al. 2011, ApJ, 743, 193 

Takami, M., Karr, J. L., Koh, H., et al. 2010, ApJ, 720, 155 

Teixeira, P.S., McCoey, C, Fich, M., & Lada, C.J. 2008, MNRAS, 384, 71 

Velusamy, T., Langer, W.D., & Marsh, KA. 2007, ApJ, 668, L159 

Werner, M.W., Roellig, T.L., Low, F.J., et al. 2004, ApJS, 154, 1 

Whitney, B. A., Wood, K., Bjorkman, J. E., & Cohen, M. 2003a, ApJ, 598, 1079 

Whitney, B. A., Wood, K., Bjorkman, J. E., & Wolff, M. J. 2003b, ApJ, 591,1049 

Wouterloot, J.G.A., & Brand, J. 1999, A&AS, 140, 177 

Wright, E. L., Eisenhardt, P. R. M., Mainzer, A. K., et al. 2010, AJ, 140, 1868 
Ybarra, J.E., & Lada, E.A. 2009, ApJ, 695, L120 (YL09) 
Ybarra, J. E., Lada, E. A., Balog, Z., et al. 2010, ApJ, 714, 469 

A. Appendix : notes on individual jets 

Below, we present a short description of both the identified jets and their candidate 
exciting sources. The numbering is that given in Tables [T] and |3j References to papers 
where possible associations are described (e.g. sub-mm cores, BLAST sources, near-IR 
imaging) are given only for Jet 1 and no longer repeated. 



This manuscript was prepared with the AAS L^Tj^X macros v5.2. 
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- Jet 1 - Extended bipolar jet composed of eight relatively bright knots. Its morphology 
is reminiscent of a S'-shape. Blue- and red-shifted jet components correspond to the 
NW (knots 1-2-3-4) and SE (knots 5-6-7-8) lobe, respectively. It is located in a region 
dominated by the strong diffuse emission from a compact HII region (263.619-0.533 

- Caswell & Haynes 1987) associated with an IR cluster (Massi et al. 2003) named 
IRS16 following the nomenclature by Liseau et al. (1992). The potential driving 
source emerges only in the IRAC bands up to 24 /im showing a SED with a Class I 
spectral index. This source is also associated with a BLAST (Olmi et al. 2009) and a 
1.2 mm dust core (Massi et al. 2007). It is not aligned with the sequence of identified 
knots, but only with few of them of the NW lobe 1,2,3,4). Therefore, we can not 
confirm this source as the jet ES. See Figure|4} 

- Jet 2 - One of the nicest examples of very compact (0.13 pc) bipolar jet with well 
defined blue- (NE) and red-shifted (SW) lobes bracketing a very young exciting 
source. This latter is point-like and detected only in the MIPS 24,70 /zm bands 
without any IRAC counterpart. Given its very red [24]- [70] color of 7.4, the system 
appears a suitable target to investigate the very early stages of the mass ejection 
(and accretion) process. An observational benefit is offered by its peculiar location 
rather isolated (see Figure[T]) and unperturbed by excessive turbulence occurring in 
proximity of HII regions or IR clusters. Remarkably, the source is associated with a 
BLAST source and with a compact dust core unresolved at 1.2 mm. See Figure[5j 

- Jet 3 - The morphology is very similar to that described for jet 2: both show 
surprisingly the same PA orthogonal to the Galactic magnetic field. It is composed 
of only two knots (the blue-shifted one at NE and red-shifted one at SW) aligned in 
a very compact (0.08 pc) and collimated shape. Both knots are symmetrical with 
respect to the position of the driving source which, again, is detected only in the MIPS 
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bands with a [24]- [70] color of 6.0. Also this jet lies in a quiet portion of VMR-D. See 
Figure[6j 

- Jet 4 - It is one of the two longest (0.68 pc) jets identified in the present survey, 
composed of the brightest detected knots. It preserves a high degree of collimation 
along all its length and shows bow-shock morphologies in both lobes oriented in 
opposite directions. Blue- and red-shifted lobes (northward and southward, composed 
by knots 1,2,3,4 and 5,6 respectively) are emanated from a well identified IR cluster 
(IRS20), where the source #98 (Massi et al. 1999) seems associated with a bipolar IR 
nebula oriented along the same jet axis. This source also coincides with a BLAST 
source and a resolved 1.2 mm dust core, with a spectral index typical of a Class I 
source. The present observations allow us to reconsider our early suggestion coming 
from a lower sensitivity survey (Lorenzetti et al. 2002b), according to which the 
northern bow shock was associated with a faint (unrelated) star. See Figure[7j 

- Jet 5 - Jet composed of well aligned knots. The membership of knot 1 to this jet is 
however doubtful, because this knot is around 170arcsec (i.e. 0.57 pc) distant from 
the others. Two objects appear as possible candidate ES : one (5a, to which we give 
higher priority) located between knots 2-6 and 7, in which case knot 1 should be 
unrelated to jet 5, and one (5b) between knot 1 and all the others, in which case the 
jet should be parsec-scale. Source 5a is barely detected with IRAC, while it is bright 
at 24 /im. Source 5b is detected at all the wavelengths and it is associated both with a 
BLAST source and with an unresolved dust core. It shows a SED with a flat spectral 
index, which possibly signals it is not a very young object. See Figure[8j 

- Jet 6 - Jet composed of five knots perfectly aligned with a candidate driving source 
north of them; this latter is classified as Class I object with a bright IRAC/MIPS 
counterpart coinciding with a 1.2 mm peak of dust emission; it is also part of the young 
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IR cluster IRS19 (Liseau et al. 1992). According to this scenario the jet morphology 
results to be one-sided with its counter-jet possibly propagating northward deep inside 
the cloud. See Figure[9j 

Jet 1 - Jet discovered in this survey. It is composed of two IRAC knots well aligned 
with a candidate exciting source with a SED with flat spectral index and having valid 
fluxes up to 70 /xm. The jet is apparently one-sided, displaying just a SE lobe. See 



Figure 10 



Jet 8 - Two knots aligned in positions symmetrically opposed with respect to the 
candidate ES located at the system center. According to our analysis the red-shifted 
(blue-) lobe extends northward (southward). The source displays a Class I SED with 
valid detected counterparts from 1.6 /im to 70 /mi; it is also associated with both a 



BLAST source and a resolved 1.2 mm dust peak. See Figure 11 



Jet 9 - Compact jet composed of five knots not perfectly aligned with the candidate 
ES. This is classified as a Class I source coinciding with a 1.2 mm dust core. See 



Figure 12 



Jet 10 - Tentatively considered as a bipolar jet emanated by a young object (whose 
SED is classified as Class I). It is associated with both a BLAST source and an 
unresolved dust peak. Such a candidate ES is located at the jet center, but slightly 
displaced with respect to the knot axis. For this reason, a different origin of the 



observed knots can not be ruled out. See Figure 13 



Jet 11 - It appears as the most compact jet of the present survey (0.07 pc), composed 
of two knots, being the southeastern the blue- and the northwestern the red-shifted 
one. They seem symmetrically positioned at both sides of a faint candidate driving 
source detected only at 24 /mi and without any association with cores at longer 
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wavelengths. Further observations are needed to better understand the true nature of 



this source. See Figure 14 



Jet 12 - Three knots very close to those of Jet 11, but presenting a different 
alignment. Blue- and red-shifted lobes are not straightforwardly recognizable, with 
some indication in favor of the blue (red)-shifted lobe displaced NE (SW). Their 
position along a well defined axis is compatible with a candidate ES located roughly 
in the central position. This source is associated with both a BLAST and 1.2 mm 



core. See Figure 15 



Jet 13 - Together with Jet 4 it appears as the longest jet of our survey. It is composed 
of six detected knots, of which knot 3 is practically coincident with the candidate ES 
while knot 4 extends slightly displaced toward SW. While knots 3-4 represent the 
most compact part of the jet, pairs of knots 1-2 and 5-6 are located far away from the 
exciting sources, but still preserving the same good alignment. In particular, both 
pairs display a shape well associated to opposite bows: the overall jet morphology 
suggests the occurrence of two (at least) episodes of matter ejection and that the 
time elapsed between them has been rather long. At a good level of confidence, we 
consider the southern (northern) lobe as the blue- (red-) shifted one. The recognized 
exciting source is bright at 24 fim, has an extremely red MIPS color ([24]- [70] = 7.9), 
and is associated with a BLAST and 1.2 mm core. Although it is a suitable candidate 
of a genuine protostellar object in a very early evolutionary stage, some doubts exist 
about its recognition as the jet ES (at least as far as the two extreme groups of knots 
is concerned). Indeed, while the source appears as one of the coolest and youngest 
of our survey, the jet (tentatively associated with it) is definitely the longest: such 
an occurrence appears rather contradictory if we assume the [24]- [70] color clS cLll ct£TG 
indicator. Indeed, if the jet were composed only of knots 3 and 4, its length should be 
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only 0.11 pc. Therefore, we prefer to not consider this jet in the analysis of Figure 22 
which is based on the estimate of the projected jet length. See Figure[T6| 



Jet 14 - Jet composed of six knots roughly aligned in E-W direction. One of them 
(indicated as 14-1) could be alternatively attributed to the Jet 15 (knot 15-5). 
Two different candidate ES have been found (labeled as 14a and 14b): the former 
(classified as a Class I, to which we give higher priority) belongs to an IR cluster and 
has a 2MASS and a MIPS counterpart, while the latter is a very faint IRAC source 



(without any WISE confirmation) with a flat spectral index. See Figure 17 



Jet 15 - We have presented this jet along with the candidate ES in previous papers 
(Lorenzetti et. al 2002; Giannini et al. 2005) to which the reader is referred to have 
a more complete view of that region. The correspondence of the knots identification 
given in those papers with that given here, is the following: the complex of knots 
previously identified as A1-A2-A3, corresponds here to 15-1, B1-B2 to 15-2, Gl to G5 
to 15-3, I to 15-4. Knot 15-5 is identified here for the first time, but its connection 
with jet 15 or jet 14 is doubtful (see also Jet 14). In the present paper other sparse 
knots (Kl to K4) have been discovered in this region, but any association with a 



plausible ES cannot be provided. See Figure 18 
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Table I. Jets found in VMR-D with IRAC. 



Jet Id. 


Number of 


length 


P.A. 


References' 1 




IRAC knots 


(")/(pc) 


(°) 




1 


8 


89/0.30 


163 


1,2 


2 


4 


34/0.13 


57 


1,2 


3 


2 


24/0.08 


57 


1 


4" 


6 


200/0.68 


174 


1,2,3 


5<= 


6(7) 


106/0.36 (321/1.09) 


48 (54) 


1,2 


6 d 


5 


74 e /0.25 e 


178 


1,2,3 


7 


2 


38/0.14 


121 




8 


2 


112/0.37 


5 


2 


9 


5 


56/0.19 


5 


2 


10 


2 


106/0.36 


43 


2 


11 


2 


20/0.07 


146 


2 


12 


3 


86/0.29 


112 


2 


13 c 


6(2) 


206/0.70 (30/0.11) 


173 


2 


14 c 


6(5) 


155/0.53 (76/0.26) 


94 (96) 




15 c ' f 


4(5) 


128/0.44 (139/0.48) 


155 (156) 


3,4 



a 1- Giannini et al. 2007, 2- Dc Luca ct al. 2007, 3- Lorcnzctti ct al. 2002, 
4- Giannini et al. 2005. 

b IRS20 in (3). 

c The number in parenthesis indicates that a different number of knots could 
belong to the jet. The jet properties are reported accordingly in the subsequent 
columns. 

d IRS19 in (3). 

c Half size than quoted in (1), see text. 
f IR.S17 in (4). 
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Table 5. Differential extinction between jet lobes. 



Jet Id. Blue knots Red knots a AAy 

(mag) 



1 


1-2-3-4 


5-6-7-8 


6 


2 


1 


4 


11 


3 


1 


2 


5 


4 


1-2-3-4 


5 


6 


Sa^ 


7 


2-3-4-5-6 


2 


5b b 


2-3-4-5-6-7 


1 


2 


8 


2 


1 


8 


!) 


5 


1-2-3 


2 


10 


2 


1 





11 


2 


1 


6 


12 


1 


2 


4 


13a c 


5-6 


1-2 


4 


13a c 


4 


3 


10 


15 


1-2 


3-4 


2 



a Knots used to compute the differential extinc- 
tion (sec Tableland Figurcs[4|to[l6). 

b For jet #5, both the possibilities that sources a 
and b are ES are considered. If the ES is source 5a, 
the pairs of symmetric knots arc 2-6 and 7 (in which 
case knot 1 does not belong to the jet), while if the 
ES is source 5b, the pairs of symmetric knots arc 1 
and 2-7. 

c For jet #13, both the possibilities that knots 
1-2/5-6 belong or not to this jet arc considered. 
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Table 6. Source parameters 01 . 



Source name 


distance 


disk incl, 


Ay (inter) 


A v (circum) 


T* 




Renv 


M onv 


L 




(pc) 


(deg) 


(mag) 


(mag) 


(K) 


(M ) 


(10 3 AU) 


(M ) 


(L©) 


1 


1000 (758) 


20 (40) 


17 (21) 


738 (1050) 


3923 (3865) 


0.87 (0.74) 


15.8 (13.2) 


43 (47) 


19 (13.7) 


2" 


500-1000 








51 








0.3-1.3 


3 b 


500-1000 








47 








0.2-0.8 


4 


692 (660) 


55 (65) 


18 (24) 


107 (876) 


4162 (4162) 


1.98 (1.98) 


2.6 (2.6) 


0.08 (0.07) 


216 (216) 


5b 


759 (832) 


20 (20) 


11 (7) 


23 (24) 


2887 (3003) 


0.15 (0.18) 


1.3 (3.0) 


0.05 (0.14) 


1.1 (1.6) 


6 


692 (912) 


50 (40) 


61 (71) 


37 (28) 


3535 (3160) 


0.40 (0.23) 


6.9 (2.2) 


0.16 (0.03) 


16 (30) 


7 


759 (501) 


30 (50) 


43 (45) 


26 (25) 


3039 (3283) 


0.20 (0.27) 


1.1° 


0.05° 


1.4 (0.8) 


8 


759 (692) 


20 (20) 


24 (21) 


22 (80) 


4370 (3949) 


0.35 (0.68) 


7.0 (9.1) 


3.5 (8.3) 


4.7 (5.6) 


9 


692 (1000) 


20 (20) 


57 (67) 


94 (2) 


3040 (3620) 


0.19 (0.42) 


3.3 (3.8) 


2.7 (3.5) 


1.3 (2.8) 


10 


955 (871) 


20 (40) 


16 (7) 


8 (29) 


2970 (3787) 


0.18 (0.65) 


3.8 (7.9) 


0.2 (0.3) 


12 (14) 


13 6 


500-1000 








43 








0.2-0.9 


14a 


955 (661) 


90 (90) 


10 (5) 


> 1000 (> 1000) 


4386 (4503) 


3.17(2.63) 


5.7 C 


0.18 c 


109 (38) 



a For each parameter, the value given in boldface refers to the best-fit model, while the value in parenthesis refers to the model whose \ 2 i s 50% higher 
than the minimum x 2 ■ 

k SED fitted with a grey-body function, by assuming a distance between 500 pc and 1000 pc. For these sources the fitted temperature is that of the 
envelope. 

c uncertain since the SED is measured up to A — 70 (im. 
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Fig. 19. — Ratio of Li_os(i) over IRAC luminosity in different bands vs. gas temperature. 
LTE approximation is adopted. 
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10 -10 

R.A. offset (arcsec) 



Fig. 20. — 12 CO(3-2) map of jet # 3 superposed on the IRAC three-color image. Contours 
are (from right to left) 1.2, 1,8, 2.4, 3.0 and 3.6 K km s _1 for the red-lobe and 0.6, 1.2, 1,8, 
2.4, 3.0 and 3.6 K km s -1 (from left to right) for the blue-lobe, respectively. The velocity is 
integrated between -3.6 km s~ l and +0.8 km s _1 for the blue-wing, and between +0.8 km 
s _1 and +3.4 km s" 1 for the right-wing. The positions of the exciting source and the H 2 
knots are evidenced as well (ellipse and boxes, respectively). 
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0.1 1 10 100 1000 0.1 1 10 100 1000 0.1 1 10 100 1000 




Fig. 21.— Fits with the Robitaille et al. (2007) model of the SED of the ES. The filled circles 
show the observed fluxes, while triangles are upper limits. BLAST data points (represented 
as upper limits) are not included in the fit. The black line shows the best-fit, and the gray 
lines show subsequent good fits. The dashed line shows the stellar photosphere corresponding 
to the central source of the best-fit model, as it would appear in the absence of circumstellar 
dust (but including interstellar extinction). 
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Fig. 22. — Projected jet length (in pc) vs. the [24]- [70] color. Arrows are upper limits. 
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Fig. 23. — IRAC line cooling vs. ES bolometric luminosity. The equation of the best-fit 
through the data points is given in the top side of the plot. Horizontal dotted lines indicate 
the range of bolometric luminosity given in Table|6j in particular the value associated with 
the best-fit model is indicated with a dot. For source # 13, the vertical line indicates Z/irac 
in case the jet is composed alternatively of two or six knots (see Appendix[A|. 



